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O : ABSTRACT 

(N 



High-resolution Hubble Space Telescope imaging observations show that the 

radial distribution of the field-decontaminated sample of 162 'blue straggler' stars 

(BSs) in the II.TIq^ Gyr-old Large Magellanic Cloud cluster Hodge 11 exhibits 

a clear bimodality. In combination with their distinct loci in color-magnitude 

Om ■ space, this offers new evidence in support of theoretical expectations that suggest 

different BS formation channels as a function of stellar density. In the cluster's 

p ■ color-magnitude diagram, the BSs in the inner 15" (roughly corresponding to the 

Q ■ cluster's core radius) are located more closely to the theoretical sequence resulting 

from stellar collisions, while those in the periphery (at radii between 85" and 100") 

ci ■ are preferentially found in the region expected to contain objects formed through 

binary mass transfer or coalescence. In addition, the objects' distribution in 

K^ color-magntiude space provides us with the rare opportunity in an extragalactic 

CN ■ environment to quantify the evolution of the cluster's coUisionally induced BS 

population and the hkely period that has elapsed since their formation epoch, 

■^ ■ which we estimate to have occurred ~4-5 Gyr ago. 

^' 

O 

cn ■ Subject headings: blue stragglers — Hertzsprung-Russell and C-M diagrams — 

stars: kinematics and dynamics — Magellanic Clouds — galaxies: star clusters: 

individual (Hodge 11) 



1. Introduction 

In globular cluster (GC) color-magnitude diagrams (CMDs), single stars that are more 
massive than the stars located at the main-sequence turn-off (MSTO) are expected to evolve 
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off the MS. If all member stars in a GC were to evolve in isolation and without under- 
going any dynamical interactions, the MS locus beyond the MSTO should be unoccupied. 
In reality, however, blue straggler stars (BSs) — which are more massive than stars at the 



MSTO — nevertheless occupy the bright MS extension (e.g., ISandagd Il953t IStrykerlll993l ). 



Two scenarios for BS formation have been suggested, involving either mass transfer be 
tween or coalescen c e of close binary 
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Lu et al.ll2010l 120111). or direct stellar collisions (e.g.. 
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The relative importance of these formation channels is still unclear. iKnigge et al.l (|2009[ ) 
found a strong correlation between the number of BSs in Galactic GC (GGC) cores and the 
associated core mass, but only a weak correlation with the predicted number of BSs formed 
through c ollisions; such a cor relation may only be seen for extremely dense GGCs, however. 
Similarly, iPiotto et al.l ( 120041 ) did not find any significant correlation between the global BS 
frequency and the collision rate in their sample clusters. However, they found that the BS 
luminosity function is noticeably different in clusters brighter and fainter than My = —8.8 
mag (for theoretical arguments, see lDavies et al.ll2004j ). Simulations indicate that the relat ive 
importance of both formation scenarios varies as a cluster evolves (JHypki fc Giersal2013l ). 



Ferraro et al.l (|2012| ) suggest that the radial profile of the relative BS frequency — the 



number of BSs normalized to that of horizontal-branch (HB) or red-giant-branch (RGB) 
stars — may be a good tracer to quantify a cluster's dynamical age. In dynamically young 
GCs, where the rela tive BS frequency is roughly fiat, e.g., u Gen, Paloma r 14, and NGG 2419 



( Ferraro et al.ll2012 ). BSs formed through mass transfer dominate (e.g., iMapelli et al.l 12006 



Hypki &: Gierszll2013l ). On the other hand, in GGs of intermediate dynarnical age the BS fre- 



quency's radial profile displays a 'bimodality' (JFerraro et al.lll997l . l2012l ; lMapelli et al.ll2006 



Hypki &: Gierszll2013l ). Stellar collisions dominate the central, crowded region, while mass 
transfer dominates the clusters' outer regions (e.g., 47 Tucanae, M3, and NGG 6752; Mapelli 
et al. 2006). In addition, BSs in dynamically old GGs (e.g., M30, M75, M 79, and M80; 



Ferraro et al. 2012) exhibit even clearer signatures of coUisional formation. iFerraro et al. 



(120091 ) discovered two distinct BS sequences in the G MD of M3 , whic h trace the theoretical 
single-age stellar-collision and mass-transfer tracks. ISills et al.l (120021 ) showed that collision 
products are not chemically homogeneous, so that coUision ally formed BSs will be redde r 
and fainter than their fully chemica lly mixed counterpa rts (JBailyn fc Pinsonneaultl Il995l ) . 
This is consistent with the results of iFerraro et al.l (120091 ). 



One may thus expect for both (dynamically) intermediate-age and old GGs that the 
number of BSs formed through either channel should be comparable. In the absence of any 
clear substructures, coUisionally induced BSs are expected to be preferentially located in a 
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GC's inner regions compared with BSs formed through mass transfer. Given the prediction 
that BSs of different origins will occupy different parts of color-magnitude space, the distri- 
bution of the different BS types in the CMD will depend on their radial distance from the 
cluster center (JFerraro et al.ll2009[ ). However, since the BSs in a given GC are not expected 
to all forr n simultaneously , the theoretical l y exp ected single-age BS sequences will be broad- 
ened (cf. iTian et al.ll2006l : ILu et al.ll2010l I2OIII ). and both BS types may partially overlap 
in color-magnitude space. However, provided that both populations are not fully mixed, a 
robust detection of two subgroups in the CMD should still be possible. 



Studies of extragalactic BSs are relatively rare. IShara et al.l (Il998l ) detected 42 BS 
candidates in the Small Magellanic Cloud cluster NGC 121 and confirmed 23 s tars as genuine 
BSs. I n the Large Magellanic Cloud (LMC) clust ers NGC 1466 a nd NG C 2257. 1 Johnson et al. 
( I1999I ) found 73 and 67 BSs, respectively, while iMackey et al.l ( l2006l )'s analysis of the LMC 
cluster ESQ 121-SC03 yielded 32 BSs. 

Using high-resolution Hubble Space Telescope (HST) observations, here we investigate 
the radial dependence of the BS frequency in the LMC GC Hodge 11. We detect a significant 
number of BSs and a clear bimodality in the cluster's radial BS frequency. We show that the 
BSs in the inner region are located more closely to the theoretical 'collision sequence' (bluer 
and fainter on average), while those in the cluster's outskirts are preferentially found in the 
'binary region' (redder and brighter). This result, combined with the discovery that Hodge 
11 hosts one of the largest BS populations in any Galactic or extragalactic cluster known to 
dateill offers unprecedented evidence in support of the theoretical predictions and allows us 
to quantify the evolution of the cluster's coUisionally formed BS population. 



Data 



Our data sets were obtained with the HSTs Wide Field and Planetary Camera-2, em- 
ploying the F555W (equi valent and hencefort h referred to as V) and F814W (~ /) filters (for 
observational details, see lde Grijs et al.ll2002[ ). The total exposure times of the observations 
used were 435 and 960 s in the V and / filters, respectively. We applied point-sp read-function 
fitting to obtain the cluster's stellar ph otometric catalog (|Dolphinll2000l . 120051 ). The cluster 
center was determined as described by Ide Grijs et al.l ( l2013l ) . We also obtained the stellar 
catalog of a nearby field region (for observational details, see lCastro et al.ll200ll ) to statisti- 



^ Although CO Cen and NGC 2419 host 313 and 266 BSs, respectively JFerraro et al.l booet 



Dalessandro et al. 2008 ). whe ther or not these objects are genuine GCs is as yet unresolved (e.g., 



Mackev fc van den Berghll2005l and references therein 
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cally reduce contamination by background stars (for procedural details, see lHu et al.ll2010l ): 
the results are shown in Fig. [1] The number of objects located in the BS region (see below 
for definition) in the original Hodge 11 CMD is 338; the decontaminated CMD contains 162 
stars in the same region. This latter sample forms the basis of our analysis in this Letter. 
Our photometric catalog is > 95% complete in the relevant magnitude range, even in the 
highest-density central region. Figure [2^ displays the spatial distribution of all observed 
stars in the Hodge 11 area covered by our HST field of view, as well as its BS population. 
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Fig. 1. — (top le ft) Original Hodge 11 CMD. (top middle) CMD of a carefully selected nearby 
field region (cf. ICastro et al.ll200ll ). (top right) Hodge 11 CMD after field-star decontam- 
ination, (bottom) Corresponding CMD zooms focusing on the BS-dominated region. The 
number of BSs and the boundaries adopted for the BS region (red outhnes; see text) are also 
indicated. The small difference between the total number of objects in the BS region (bot- 
tom left) and the sum of the BSs identified in the bottom-middle and botto m-right panels i s 
due to the statistical nature of our field-star decontamination procedure (cf. IHu et al.ll2010l ). 



We divided the MS and RGB regions in the Hodge 11 CMD into 32 bins of ~ 0.30 mag 
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Fig. 2. — (left) Spatial distribution of stars observed in the Hodge 11 area. Red point: 
Cluster center. Blue points: BSs. Red (partial) concentric circles indica te radii of -R = 15, 85 , 
and 100". (right) Corresponding HST CMD, including the best-fitting iBertelli et all ^M) 
isochrone (blue). Red points represent our BS selection; the inset shows the photometric 
uncertainties for all BS candidates. Black points with error bars have photometric errors 
AV > 0.1 mag (color uncertainties have been omitted for clarity). Thick solid lines delineate 
the BS region adopted. 



each between V = 18.3 and 27.0 mag. In each magnitude range, we determined the peak and 
3a spread of the color distribution using a Gaussian function (which provided appropriate fits 
at any magnitude). Spectroscopic abundance deterrn i nations show that Hod ge 11 is metal- 
poor, [Fe/H] = -2.1 ± 0.2 dex jCowlev et all Il982l : loiszewski et~ZI Il99l[ ) to -1.86 dex 
JGrocholski et al.lboOBJ l. Adopting a metaUicity Z = 0.0002 ([Fe/H] ^ -1.98 dex) and using 



the ridgeline thus obtained, we determined the best-fitting isochrone based on the Padova 
stellar evolutionary models (JBressan et al.l 120121 ) and calculated the M STO locus. We thus 



determined a cluster age of log(t yr ) = 10.07 ±0.01 (11.7_o;i Gyr; cf. iJohnson et al.lll999l ) 



and an e xtinc tion of E{B — V) = 0.09 ± 0.01 mag. The latter is in excellent agreement with 
Walked (jl993l )'s determination of E{B — V^) = 0.08 ± 0.02 mag. We adopted the canonical 
LMC distance modulus, (m — M)q = 18.50 mag, so that 1" = 0.24 pc. 

Many of the cluster's HB stars are located in the extreme blue HB (EHB). Our sample 
BSs can easily be distinguished from these (E)HB stars by adopting a minimum (E)HB-BS 
difference in color and magnitude of 0.3 mag. Our selection is robust, since the (E)HB stars 
and BSs are clearly separated in the CMD (see Fig. [2]) and the absolute photometric errors 
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are <C 0.3 mag. 

Finally, we identify a star as a likely BS if it meets all of the following criteria: 

1. It is bluer and brighter than the single-star MSTO locus. 

2. Its distance in color-magnitude space from the MSTO and MS ridgeline is > 3a {a 
refers to the intrinsic spread in MS color at the magnitude of interest). 

3. Its photometric uncertainty is < 0.1 mag in both filters. 

Figure[2]D shows the Hodge 11 HST CMD, where the inset focuses on the region typically 
populated by BSs. The red points represent our (field-decontaminated) BS candidates. The 
small black points with error bars are objects whose photometric errors are > 0.1 mag, while 
most of the error bars of the red points are smaller than the symbol size. Inspection on the 
original images of the sources with large error bars revealed that they either are affected 
by instrumental artefacts (including low signal-to-noise ratios or closely juxtaposed to bad 
pixels) or may be unresolved blends. We thus discarded these objects from our BS sample. 



3. Analysis and results 

Figures [3fe' and [3b show the cluster's stellar number density as a function of radius for 
all stars with V < 2'b mag. Based on the catalog of the nearby field region, we calculated 
the field's average density value. We adopted the distance where the cluster's monotonically 
decreasing number- density profile reaches the average field level, R = IO2.3II3 arcsec, as 
the cluster region. 

In addition to our BS sample, we selected a sample containing RGB stars covering the 
magnitude range from V = 19.60 to 22.05 mag. We normalized the number of BSs to the 
number of RGB stars: /bss = A^bss/A^rgb, where A^bss and ./Vrgr represeri t the number of 



BSs and RGB stars, respectively. Previous studies (e.g.. lFerraro et al.ll2003l ) normalized the 
number of BSs to either or both HB or/and RGB stars. We adopted RGB stars as reference 
for reasons of consistency with previous analyses. We subsequently determined the radial 
dependence of the BS fraction, /bss(-R): see Fig. ^. Normalization of the BS number to the 
number of HB stars results in a similar, although somewhat flatter trend. This BS-fraction 
profile is very simila r to the equivalent radial dependences in 47 Tuc, M3, and NGC 6752 



flMapeUi et al.ll2006h . 



The number of BSs in Hodge 11 is extremely large. iFerraro et al.l (120071 ) summarized 



the number of BSs detected in 26 GGCs, ranging from four (NGC 5024) to 137 (NGC 5272). 
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In Hodge 11, we detected 150 BSs within R = 100"; the total number of BSs on all four 
chips of our Hodge 11 observations reaches 162 (see Section H] for an assessment of the level of 
completeness of our observations' areal coverage). Since the BS distribution associated with 
Hodge 11 is very extended, the number of BSs located in our outer subsample, R G [85", 100"] , 
is significant. This annulus contains 27 BSs, and hence this large number of BSs strengthens 
the statistical reliability of our conclusions. 



o 4 

CD ^ 
CO 

03 2 



0.14 

CD 

1 0.12 

^. 0.1 

c?0.08 
10 



50 100 

R (arcsec) 



1 



u 

CD 
CO 



bJD. 

O 



15 



"^ 



1.5 

logi? [arcsec] 



I I I I I 

20 25 30 35 

R (arcsec) 



I I I 



40 



45 



50 



55 



Fig. 3. — (a) Hodge 11 number- density profile (corrected for partial areal coverage). The red 
dashed line indicates the average field density, implying a cluster size oi R = 102.3l2'3 arcsec. 
(b) Corresponding logarithmic number-density profile, (c) BS-fraction profile (normalized 
to the number of RGB stars) of Hodge 11, including Poissonian uncertainties. 



Mackey &: Gilmord (120031 ) determined a cluster core radius of -Rcore = 12.14 ± 0.66 
arcsec. Our number-density profile exhibits a relatively sharp radial decline to R ~ 15- 
20". Therefore, we defined an inner BS subsample covering the cluster core at R < 15" 
(Fig. m left). The red and blue points in the CMD of Fig. H] (left) represent the inner 
{R < 15") and outer {R G [85", 100"]) subsamples. Both contain 27 BSs. A distinct difference 
is seen between the mean loci of both samples in color-magnitude sp ace. We also show 



the separation between the two BS subpopulations in M30 found by (JFerraro et al.l 12009 



dashed hne). These authors suggested that BSs below this 'critical' line might preferentially 
result from stellar collisions in relatively high-density regions, while BSs above the line are 
predominantly thought to have formed through mass transfer between binary companions 
in lower- density areas. 



Considering that M30 and Hodge 11 may have different evolutionary histories and en- 



vironments, using the M30 critical line to analyze the color-magnitude distribution of the 
Hodge 11 BSs may not be appropriate. If we shift the zero-age main sequence (ZAMS) up- 
ward in the CMD by 0.75 mag — which represents the ZAMS of equal-mass binary systems 
(B-ZAMS) — we see that this provides an excellent se paration between b oth subsamples. 
Hence, our results for Hodge 11 support the claims of iFerraro et al.l (120091 ) for M30: Only 
one outer-sample BS is unambiguously located in the 'bottom region' (labeled 'A' in Fig. HI 
middle), while the 26 remaining outer-sample BSs are all located above the B-ZAMS (mod- 
ulo the photometric uncertainties). In contrast, for our R < 15" inner-sample BS selection, 
only four BSs are unambiguously located in this CMD region, given the photometric errors. 
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Fig. 4. — (left) Color-magnitude distribution of BSs in Hodge 11. Red points: BSs at 
R < 15"; blue points: BSs at i? G [85", 100"]. The black dashed line indicate s approximately 
the critical line separating the two BS sequences in M30 (IFerraro et al.ll2009h : the black solid 
line indicates the predicted ZAMS adjusted by —0.75 mag. (right) Evolutionary tracks of 
stars with masses from 1.0 to 1.6Mq in steps of O.IM© (bottom to top; Bertelli et al. 2008) 
for Z = 0.0004 (closest available metallicity). Different colors indicate the evolutionary 
timescales for stars to evolve from the ZAMS to their current positions, in units of 10^ yr. 
Red open circles: Inner-sample BSs {R < 15") located above the B-ZAMS. 



Discussion 



Since the cluster's central region is very crowded, the BSs in this region are most likely 
preferentially the result of stellar collisions. On the other hand, since the low number den- 
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sity at large radii renders a high frequency of direct stellar collisions unlikely, the BSs in the 
cluster's periphery most likely originate from mass transfer between the individual members 
in binary systems. Based on the B-ZAMS adopted, the mixture of the two BSs subsamples 
in the CMD seems 'unbalanced:' only few outer-sample BSs are found in the bottom re- 
gion, while inner-sample BSs frequently 'pollute' the top region. We suggest three possible 
explanations for this observation: 

1. The spatial distribution of the BSs in Hodge 11 is a two-dimensional projection of a 
3D spatial distribution. Some apparent inner-sample BSs may not be genuine inner- 
sample members, but could be outer-sample BSs that are projected close to the cluster 
center. On the other hand, outer-sample BSs must always be located truly far away 
from the center. Geometry considerations, assuming a uniform BS distribution in the 
outer annulus, R G [85", 100"], suggest that 4±2 (Poissonian uncertainty) outer-sample 
BS may be projected along the full line of sight through the cluster onto radii R < 15". 

2. BSs evolve and trace paths in the CMD. Without any new djTiamical interactions, the 
BSs in the bottom region will move upward, but the BSs in the top area will never 
move downward. To estimate the approximate time necessary for BSs to evolve from 
the ZAMS to the region above the B-ZAMS, in Fig. H] (right) we show evolutionary 



tracks from iBertelli et al.l (120081 ) for stars of different masses. Below, we select those 
inner-sample BSs that are located above the B-ZAMS (red open circles) and compare 
their loci to the tracks for stars with masses from 1.0 to I.QMq (bottom to top; I.QMq 
is approximately 2Mmsto)- Different colors indicate different evolutionary timescales 
for stars to evolve from the ZAMS to their current positions. 

3. There may still be some primordial binaries in the inner region, which may merge to 
form BSs. Once formed, they will be located in the top region. 

One outer-sample BS, object 'A' in Fig. H] (middle), is found in the CMD's bottom 
region. Its locus cannot have been caused by the expected evolutionary color spread. This 
star could either be a projected foreground star or a collisionally formed cluster BS tha t 



may have been ejected dynamically to the cluster's periphery (cf. ISigurdsson et al.lll994l ). 
Such ejected BSs will subsequently sink to the center again through the effects of two-body 
relaxation (dynamical mass segregation). Compared with the other cluster BSs, this star is 
relatively faint (and hence of relatively low mass), so that this relaxation process is likely 
still underway. 

Accounting for scenarios 1-3 above, the radial distance bias of BSs in the Hodge 11 CMD 
is nevertheless significant, which suggests that our result is robust. In addition, compared 
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with M30, the number of field BSs detected in the Hodge 11 field at radii beyond the 
cluster's outer radius is very large. We identified 12 BSs in the partially observed annulus 
R G [100", 130"], whose ffST-based coverage is only 9.3% complete (for comparison, our 
coverage of the actual cluster area, R < 100", is 50.2% complete). Nevertheless, if we 
adopt an outer boundary to the outer-sample BS population such that all observed BSs are 
included and compare their color-magnitude distribution with that of the inner-sample BSs, 
the distinction between the two samples remains significant. 

If we treat all selected objects as genuine cluster BSs, the locus of the most massive BS — 
which should originally have been located in the 'bottom region' of the CMD but is found 
above the B-ZAMS — may have been caused by stellar evolution. Geometry considerations 
imply that 4 ± 2 stars may be seen in projection, although nine inner-sample BSs that may 
have crossed the B-ZAMS are actually detected. This means that 5 ± 2 objects could be 
evolved BSs. If correct, this allows us to estimate their time of birth at 4-5 Gyr ago. The 
BSs characterized by magnitudes between V = 21.0 and 21.5 mag are probably also evolved 
stars. For BSs fainter than K = 21.5 mag, the time required to evolve to their current loci 
is so long that their positions in the CMD may be owing to 2D projection of the cluster's 
3D spatial distribution. 



5. Conclusion 

Based on the predictions of stellar evolutionary models, BSs resulting from direct stellar 
collisions and from mass transfer between binary companions will be located along two 
distinct sequences. However, two well-separated sequences can only be observed if all BS 
formation in a given cluster occurred recently or if all BSs are short-lived. Long-term, 
continuous BS formation will render the expected boundary less distinct. Nevertheless, it 
is expected that BSs formed through either of these two mechanisms show the signature of 
a clustercentric radial dependence in their color-magnitude distribution. We analyzed the 
HST CMD of the GC Hodge 11, where we found a clear signature of such a spatial bias. The 
average BS in the inner region is relatively faint, while BSs in the outer region are brighter. 
Even when considering the possible effects of projection and evolution, this difference is 
robust. Our result hence offers strong evidence in support of the theoretical expectation of 
dual-mode BS formation and allowed us to quantify their evolutionary timescales. 
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